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Abstract. This work is the continuation of the recent paper [9] devoted to the density- 
dependent incompressible Euler equations. Here we concentrate on the well-posedness issue 
^— ■>) , in Besov spaces of type B^^^ embedded in the set of Lipschitz continuous functions, a func- 

■ tional framework which contains the particular case of Holder spaces and of the endpoint Besov 

space B^a^i. For such data and under the nonvacuum assumption, we establish the local well- 
posedness and a continuation criterion in the spirit of that of Beale, Kato and Majda in ^ . 

In the last part of the paper, we give lower bounds for the lifespan of a solution. In dimension 
two, we point out that the lifespan tends to infinity when the initial density tends to be a 
constant. This is, to our knowledge, the first result of this kind for the density-dependent 
incompressible Euler equations. 



< 

^ I 1. Introduction and main results 

This work is the continuation of a recent paper by the first author (see ^) devoted to the 
density- dependent incompressible Euler equations: 

dtp + u-Vp = 0, 

p{dtu + u ■ Vn) Vn = p/, 

divti = 0. 



(N 

in 
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Recall that the above equations describe the evolution of the density p = p{t, x) G and 
of the velocity field u = u{t, x) G of a nonhomogeneous inviscid incompressible fluid. The 
O ■ time dependent vector field / stands for a given body force and the gradient of the pressure 

I vn is the Lagrangian multiplier associated to the divergence free constraint over the velocity. 

We assume that the space variable x belongs to the whole M'^ with N >2. 

There is an important literature devoted to the standard incompressible Euler equations, that 
^ , is to the case where the initial density is a positive constant, an assumption which is preserved 

during the evolution. In contrast, not so many works have been devoted to the study of ([1]) in 
the nonconstant density case. In the situation where the equations are considered in a suitably 
smooth bounded domain of or M^, the local well-posedness issue has been investigated by 
H. Beirao da Veiga and A. Valli in [3l HI [5] for data with high enough Holder regularity. In [8], 
we have proved well-posedness in with s > 1 -|- N/2 and have studied the inviscid limit in 

-+i 

this framework. Data in the limit Besov space were also considered. 

As for the standard incompressible Euler equations, any functional space embedded in the set 
C^'^ of bounded globally Lipschitz functions is a candidate for the study of the well-posedness 
issue. This stems from the fact that System ([1]) is a coupling between transport equations. 
Hence preserving the initial regularity requires the velocity field to be at least locally Lipschitz 
with respect to the space variable. As a matter of fact, the classical Euler equations have been 
shown to be well posed in any Besov space B^ ,. embedded in C°'^ (see [H HSl IB] and the 
references therein), a property which holds if and only if {s,p,r) G M x [l,-|-oo]^ satisfies 

N N 

(C) s > 1 -\ or s = l-| and r = 1 . 

p p 



Date: May 7, 2013. 



1 



2 



R. DANCHIN AND F. FANELLI 



In [9] , we extended the results of the homogeneous case to ([T]) (see also [10] for a similar study 
in the periodic framework). Under condition (C) with 1 < p < oo we established the local well- 
posedness for any data {po,uo) in Bp ,, such that po is bounded away from zero. However, we 
have been unable to treat the limit case p = oo unless the initial density is a small perturbation 
of a constant density state, a technical artifact due to the method we used to handle the pressure 
term. 

In fact, in contrast to the classical Euler equations, computing the gradient of the pressure 
involves an elliptic equation with nonconstant coefficients, namely 

(2) div (aVn) = div F with F := div {f - u ■ Vu) and a := 

Getting appropriate a priori estimates given that we expect the function p to have exactly the 
same regularity as VII is the main difficulty. In the framework and, more generally, in the 
Sobolev framework H^, this may be achieved by means of a classical energy method. This is 
also quite straightforward in the Bp^^ framework if a is a small perturbation of some positive 
constant function a, for the above equation may be rewritten 

aAU = div F + div ((a - a)Vn) . 

Now, if a — a is small enough then one may take advantage of regularity results for the Laplace 
operator in order to "absorb" the last term. 

If 1 < p < oo and a is bounded away from zero then it turns out that combining energy 
arguments similar to those of the case and a harmonic analysis lemma allows to handle the 
elliptic equation This is the approach that we used in [9]. However it fails for the limit 
cases p = 1 and p = oo. 

In the present work, we propose another method for proving a priori estimates for ([2]). In 
addition to being simpler, this will enable us to treat all the cases p G [l,oo] indistinctly 
whenever the density is bounded away from zero. Our approach relies on the fact that the 
pressure H satisfies (here we take / = to simplify) 

(3) AH = -pdiv(ii- Vii) + Vlog/9-VH. 

Obviously, the last term is of lower order. In addition, the classical theory ensures that 

a*||Vn||^2 < Vu||j;^2 with a* := inf a(x). 

Therefore interpolating between the high regularity estimates for the Laplace operator and the 
estimate allows to absorb the last term in the right-hand side of ([3]). 

In the rest of the paper, we focus on the case p = oo as it is the only definitely new one and 
as it covers both Holder spaces with exponent greater than 1 and the limit space B^^ ^ which 
is the largest one in which one may expect to get well-posedness. 

Before going further into the description of our results, let us introduce a few notation. 

• Throughout the paper, C stands for a harmless "constant" the meaning of which depends 
on the context. 

• If a = (a^,a^) and b = (6^,6^) then we denote a A 6 := a^b"^ — a^b^ . 

• The vorticity Q, associated to a vector field u over is the matrix valued function 
with entries 

fijj := dju!' — diu' . 

If iV = 2 then the vorticity may be identified with the scalar function uj := diu^ — d2U^ 
and if N = 3, with the vector field V x u. 

• For all Banach space X and interval / of R, we denote by 6(1; X) (resp. Ch{I;X)) 
the set of continuous (resp. continuous bounded) functions on / with values in X. If X 
has predual X* then we denote by Ci„(/;X) the set of bounded measurable functions 
f : I ^ X such that for any 4> G X* , the function t i— )• {f{t),(j))xxX* is continuous 
over /. 
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• For p G [1,00], the notation LP{I;X) stands for the set of measurable functions on / 
with values in X such that t 1— t- belongs to In the case / = [0,T] we 
alternately use the notation Lj,{X). 

• We denote by L^^^{I) the set of those functions defined on / and valued in X which, 
restricted to any compact subset J of /, are in L/'{J). 

• Finally, for any real valued function a over R^, we denote 

a* := inf a(x) and a* := sup a{x). 

Let us now state our main well-posedness result in the case of a finite energy initial velocity 
field. 

Theorem 1. Let r be in [l,oo] and s G M satisfy s > 1 if r ^ 1 and s > 1 if r = 1. Let pQ 
be a positive function in ^ bounded away from 0, and uq be a divergence- free vector field 
with coefficients in ^ Ci L"^ . Finally, suppose that the external force f has coefficients in 
L^{[—Tq,Tq];B^j.) n C([— Tq, Tq]; L^) for some positive time Tq. 

Then there exists a time T £]0,Tq] such that System ([T]) with initial data {pq,uq) has a 
unique solution (/?, u, VII) on [—T,T] xM^, with: 

• p in C{[—T,T]]B^j.) and bounded away from 0, 

• u in C([-r,r];S^'jnCi([-r,T];L2) and 

. vn zn Li([-r,r];i?^^jnc([-r,r];L2). 

Lf r = 00 then one has only weak continuity in time with values in the Besov space B^ ^ . 

In the above functional framework, one may state a continuation criterion for the solution to 
([1]) similar to that of Theorem 2 of [9] : 

Theorem 2. Let {p,u,VIl) be a solution to System ([1]) on [0, r*[xM^, with the properties 
described in Theorem{^for all T <T* ; suppose also that we have 

(4) / (^||Vu||loo + ||Vn||^s-i^ < 00 . 

J 

Lf T* is finite then (p, u, VII) can be continued beyond T* into a solution of ([T]) with the 
same regularity. Moreover, if s > 1 then one may replace in ([4]) the term \\Vu\\l°° with ||$7||ioo. 
A similar result holds for negative times. 

From this result, as our assumption on (r, s) implies that B^j. ^ , keeping in mind that 
-^00 1 largest Besov space included in C^'^ , we immediately get the following: 



Corollary 1. The lifespan of a solution in B^ ^ with s > 1 is the same as the lifespan in B\ 



00,1 ■ 



As pointed out in [9j, hypothesis uq G L^ is somewhat restrictive in dimension = 2 as 
if, say, the initial vorticity wq is in L^ then it implies that ujq has average over M'^ . In 
particular, assuming that uq G L'^{E?') precludes our considering general data with initially 
compactly supported nonnegative vorticity (e.g. vortex patches as in [7], Chapter 5). 

The following statement aims at considering initial data with infinite energy. For simplicity, 
we suppose the external force to be 0. 

Theorem 3. Let {s,r) be as in TheoremUl Let po G B^^ be bounded away from 0, and 
uoGB^^^.nW'^. 

Then there exist a positive time T and a unique solution (/?, u, VII) on [—T, T] x of 
System ([1]) with external force / = 0, satisfying the following properties: 

• /9 G C{[—T,T];B^j.) bounded away from 0, 

• u G C([-r,r];5^'^^nVFi'4) and dtu£C{[-T,T];L'^), 
. vn G LH[-T,T];B^^^,)nC{[-T,T];L^). 

As above, the continuity in time with values in B^ ^ is only weak if r = 00 . 
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Remark 1. Under the above hypothesis, a continuation criterion in the spirit of Theorem\^ 
may be proved. The details are left to the reader. 

Let us also point out that there is some freedom over the W^''^ assumption (see Remark^ 
below). 

On the one hand, the existence results that we stated so far are local in time even in the 
two-dimensional case. On the other hand, it is well known that the classical two-dimensional 
incompressible Euler equations are globally well-posed, a result that goes back to the pioneering 
work by V. Wolibner in jl6] (see also [17^ ITTI IT5] for global results in the case of less regular 
data). In the homogeneous case, the global existence stems from the fact that the vorticity cj 
is transported by the flow associated to the solution: we have 

dtuj + u ■ Vuj = . 

In the nonhomogeneous context this relation is not true any longer: we have instead 

(5) dtu + u-Vuj + V A Vn = . 

If the classical homogeneous case has been deeply studied, to our knowledge there is no 
literature about the time of existence of solutions for the density-dependent incompressible 
Euler system. In the last section of this paper, we establish lower bounds for the lifespan of a 
solution of 

Roughly, we show that in any space dimension, if the initial velocity is of order e (e small 
enough) , without any restriction on the density of the fluid then the lifespan is at least of order 
(see the exact statement in Theorem (H). 

Next, taking advantage of Equality ^ and of an estimate for the transport equation that 
has been established recently by M. Vishik in [15] (and generalized by T. Hmidi and S. Keraani 
in [E]), we show that the lifespan of the solution tends to infinity if po — 1 goes to 0. More 
precisely, Theorem [5] states that if 

IIpo - = and ||wo||bo + ll^^ollLa = e 

oc.l oo,l 

with r] small enough, then the lifespan is at least of order e"-*^ log(log r/^^). 

The paper is organized as follows. In the next section, we introduce the tools needed for 
proving our results: the Littlewood-Paley decomposition, the definition of the nonhomogeneous 
Besov spaces i?*^ and the paradifferential calculus, and finally some classical results about 
transport equations in Bp ^. and elliptic equations. Sections [3] and H] are devoted to the proof of 
our local existence statements first in the finite energy case and next if the initial velocity is in 
VF^'^. Finally, in the last section we state and prove results about the lifespan of a solution of 
our system, focusing on the particular case of space dimension N = 2. 

2. Tools 

Our results mostly rely on Fourier analysis methods based on a nonhomogeneous dyadic parti- 
tion of unity with respect to the Fourier variable, the so-called Littlewood-Paley decomposition. 
Unless otherwise specified, all the results which are presented in this section are proved in [T]. 

In order to define a Littlewood-Paley decomposition, fix a smooth radial function x supported 
in (say) the ball .6(0, |), equals to 1 in a neighborhood of i?(0, |) and such that r i— )• xi''' ^r) 
is nonincr easing over R_|_, and set ¥5(0 = ^(l) ~ x(C)- 

The dyadic blocks {lS.j)j^i are defined b}Q 

Aj := if j < -2, A_i := x{D) and A^ := (f{2-W) if j > 0. 



Throughout we agree that f{D) stands for the pseudo-differential operator T ^{fJ-u). 
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We also introduce the following low frequency cut-off: 

Sju:=x{2-'D)= Aj, for j > 0. 

i'<i-i 

The following classical properties will be used freely throughout in the paper: 

• for any u € S', the equality u = AjU holds true in 5'; 

• for all u and v in S', the sequence {Sj-iu Ajv)j^^ is spectrally supported in dyadic 
annuli. 

One can now define what a Besov space ^ is: 

Definition 1. Let u he a tempered distribution, s a real number, and 1 <p,r< oo. We set 



miB" 



2''-^'^|| Ajtt||^p ] if r < oo and ||n||_Bs^ := sup (2-'''|| Aju||lp) . 



We then define the space B^^ as the subset of distributions u £ S' such that \\u\\b^^. is finite. 

From the above definition, it is easy to show that for all s G M, the Besov space -B| 2 
coincides with the nonhomogeneous Sobolev space H^. Let us also point out that for any k £N 
and p G [1, 00], we have the following chain of continuous embedding: 

where W^'^ denotes the set of L'p functions with derivatives up to order k in L^. 

The Besov spaces have many nice properties which will be recalled throughout the paper 
whenever they are needed. For the time being, let us just recall that if Condition (C) holds true 
then Bp ^ is an algebra continuously embedded in the set C'^'^ of bounded Lipschitz functions 
(see e.g. [1], Chap. 2), and that the gradient operator maps B^ ^. in Bp~^. 

The following result will be also needed: 

Proposition 1. Let F : — )• M 6e a smooth homogeneous function of degree m away from a 
neighborhood of the origin. Then for all {p,r) G [l,oo]^ and s G M, Operator F{D) maps Bp ,^ 
in B'p-"^. 

Remark 2. Let V be the Leray projector over divergence free vector fields and Q := Id — V. 
Recall that in Fourier variables, we have for all vector field u 

Therefore, both (Id — A_i)'P and (Id — A_i)Q satisfy the assumptions of the above proposition 
with m = hence are self-map on Bp ^ for any s G M and 1 < p, r < 00. 

The following lemma (referred in what follows as Bernstein's inequalities) describes the way 
derivatives act on spectrally localized functions. 

Lemma 1. Let < r < R. A constant C exists so that, for any nonnegative integer k, any 
couple {p, q) in [1, 00]^ with q > p > 1 and any function u of L^ , we have for all A > 0, 

SuppSc S(0,Ai?) =^ \W^u\\li < C*-'+^A^'+^^^"9)||u||lp; 

SuppSc U e K^/rA < 1^1 < R\} =^ C-^-^\^\\u\\lv < \\V^u\\lp < C^+^\^\\u\\lv. 

The first Bernstein inequality entails the following embedding result: 

Corollary 2. The space Bp\^^_^ is continuously embedded in the space Bp'^ ,.^ whenever 1 < pi < 
P2 < 00 and 

S2 < si — N/pi + N/p2 or S2 = si — N/pi + N/p2 and 1 < ri < r2 < 00. 
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Let us now introduce the paraproduct operator and recall a few nonlinear estimates in Besov 
spaces. Constructing the paraproduct operator relies on the observation that, formally, any 
product of two tempered distributions u and v, may be decomposed into 

(6) uv = TuV + TyU + R{u, v) 

with 

TuV := ""^^ Sj^iuAjV, TyU := ""^^ Sj^ivAjU and i?(u,t;):=^^ AjuAj/v. 
i 3 j |j'-il<i 

The above operator T is called "paraproduct" whereas R is called "remainder". 

The paraproduct and remainder operators have many nice continuity properties. The follow- 
ing ones will be of constant use in this paper (see the proof in e.g. [Ij, Chap. 2): 

Proposition 2. For any (s,p,r) G M x [1,cxd]^ and t < 0, the paraproduct operator T maps 
L°° X B^ .^. in Bp J., and B^^^ x Bp ^. in -B^^*. Moreover, the following estimates hold: 

\\Tuv\\b^ ^ C\\u\\l°°\\Vv\\ ns-1 and llTuvllns+t < ClluWot ||Vw||ns-i- 

M ^ ll-^p,r — II 11-^ II ''-Dp^r " "-Dp,r — " "-'-'00,00 H "-Dp,r 

For any {si,pi,ri) and {s2,P2,r2) m M x [1, oo]^ such that si + S2 > 0, 1/p := l/pi + l/p2 < 1 
and 1/r := 1/ri + l/r2 < 1 the remainder operator R maps Bp^ j._^ x Bp^ ^.^ in Bp^^f^^. 

Combining the above proposition with Bony's decomposition ([6]), we easily get the following 
"tame estimate": 

Corollary 3. Let a be a bounded function such that Va G Bp~^ for some s > and {p, r) G 
[l,oo]^. Then for any b G Bp^. n L°° we have ab G Bp j. n L°° and there exists a constant C 
depending only on N, p and s such that 

\\ab\\Bi^ < c(^\\a\\Loo\\b\\Bs^^ + ||b||L- 1| Va||^.-i 

The following result pertaining to the composition of functions in Besov spaces will be needed 
for estimating the reciprocal of the density (see the proof in [9]). 

Proposition 3. Let L be an open interval of M and F : / — )• M, a smooth function. Then for 
all compact subset J C L, s > and {p,r) G [l,oo]^ there exists a constant C such that for all 
function a valued in J and with gradient in Bp~.^, we have V(F(a)) G Bp~^ and 

||V(F(a))||^.-x<C||Va||^.-i. 

Our results concerning Equations ([T]) rely strongly on a priori estimates in Besov spaces for 
the transport equation 



(T) 



j dta + v ■\7a = f, 

We shall often use the following result, the proof of which may be found in e.g. [1]. 



Proposition 4. Let 1 < r < oo and a > (a > —1 if divv = 0). Let oq G B^^., f G 
L^{[0,T]; B^ J.) and v be a time dependent vector field in Cb{[0,T] x M^) such that 

Vv G L\[0,T];L°°) if a < 1, 

Vv G L^{[0,T];B^J:) if a > 1, or = r = 1. 

Then Equation (T) has a unique solution a in 
• the space C{[0,T]; B^ j.) if r < oo, 

. the space (f].,^,C{[0,T]- B^)) f]CU[0,n B-) if r = ^. 
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Moreover, for all t € [0,T], we have 

(7) e-^^W||a(t)||Bs,, < llaolbg. + / e-^^(*')||/(Obg. dt' 



with V'{t) :-- 



\\Vv{t)\\L^ if a<l, 
\\Vv(t)\\^cr-i if a > 1, or a = r = 1. 
If a = V then, for all a > (cr > —1 if divv = 0), Estimate ([7]) holds with V'{t) := ||Va(t)||L>^. 
Finally, we shall make an extensive use of energy estimates for the following elliptic equation: 

(8) - div (aVn) = div F in 

where a = a{x) is a given suitably smooth bounded function satisfying 

(9) a* := inf a{x) > 0. 

We shall use the following result based on Lax-Milgram's theorem (see the proof in e.g. [9]). 

Lemma 2. For all vector field F with coefficients in , there exists a tempered distribution 11, 
unique up to constant functions, such that VII G and Equation ([8|) is satisfied. In addition, 
we have 

(10) a*||Vn||i2 < \\F\\l2. 

3. Proof of Theorem [T] 

Obviously, one may extend the force term for any time so that it is not restrictive to assume 
that To = +00. Owing to time reversibility of System ([TJ, we can restrict ourselves to the 
problem of evolution for positive times only. For convenience we will assume r < oo; for 
treating the case r = oo, it is enough to replace the strong topology by the weak topology, 
whenever regularity up to index s is involved. 

We will not work on System ([1]) directly, but rather on 

' dta + u-Va = Q 

(11) \ dtu + u-Vu + aVn = / 

- div (aVn) = div (u • VVu - f) , 

where we have set a := 1/p. 

The equivalence between ([T]) and (|lip is given in the following statement (see ^). 

Lemma 3. Let u be a vector field with coefficients in C^{[0,T] x M^) and such that Qu G 
C^([0, T]; L^) . Suppose also that VII G C([0,T];L^). Finally, let p be a continuous function on 
[0, T] X such that 

(12) 0<p* </)</)*. 

Let a := 1/p. If div u(0,-) = in then (p, it, VII) is a solution to ^ if and only if 
(a, n, Vn) is a solution to (llip . 

We now come to the plan of this section. First of all, we shall prove a priori estimates for 
suitably smooth solutions of ([1]) or (jlip . Even though those estimates are not needed for proving 
Theorem [H they will be most helpful to get the existence. As a matter of fact, the construction 
of solutions which will be proposed in the next subsection amounts to solving inductively a 
sequence of linear equations. The estimates for those approximate solutions turn out to be the 
same as those for the true solutions. In the last two subsections, we shall concentrate on the 
proof of the uniqueness part of Theorem [1] and of the continuation criterion stated in Theorem 
[2] (up to the endpoint case s = r = 1 which will be studied in the next section). 
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3.1. A priori estimates. Let (a,u, VII) be a suitably smooth solution of System pT]) with 
the required regularity properties. In this subsection, we show that on a suitably small time 
interval (the length of which depends only on the norms of the data), the norm of (a, u, VII) 
may be bounded in terms of the data. 

Recall that according to Proposition [3] the quantities ||a||B|^^, and are equivalent 

under hypothesis (jl2p . This fact will be used repeatedly in what follows. 

3.1.1. Estimates for the density and the velocity field. Let us assume for a while that divn = 0. 
Then (p, u, VH) satisfies System ([T]) and the following energy equality holds true: 



(13) 



\\Vp(t)Ht)\\l, = \\^uo\\l2 + 2 f (I pf-udx)dT. 

Jo \Jm J 



Moreover, from the equation satisfied by the density, we have that p{t,x) = pQ (^tp^^ (x)) , 
where ip is the flow associated with u; so, p satisfies p^ . Hence, from relation (fT3|) . we obtain 
the control of the norm of the velocity field: for all t G [0, Tq] , we have, for some constant 
C depending only on and p* , 

(14) \\u{t)h2 < c(\\uo\\l^ + f \\f{T)\\L^dT\ . 



\ JO J 

Next, in the general case where divu need not be 0, applying Proposition|l]yields the following 
estimates: 

(15) \\a{t)\\Bs^^^ < \\ao\\B^^_,. exp(^C \\u\\B^^ ,.dT^ 

(16) ||n(t)||ij^_^ < exp(^C \\u\\B^ ^dT^ ■ (^\uo\\b^^ ^ + 

where, in the last line, we have used the fact that ^ , under our hypothesis, is an algebra. 

Remark 3. Of course, as p and a verify the same equations, they satisfy the same estimates. 

3.1.2. Estimates for the pressure term. Let us use the low frequency localization operator A_i 
to separate VH into low and high frequencies. We get 

iivniiB- < ||A_iVn||B. + ||(id- A_i)vn||Bs . 

II 1 1 -^oo ,r — II ^ ii-^oo,r ^ ^ ' ii-^oo,r 

Observe that (Id — A_i)Vn may be computed from AH by means of a homogeneous mul- 
tiplier of degree —1 in the sense of Proposition [TJ Hence 

(17) ||(id-A_i)vn||Bg,„ < c||An||^^-i. 

For the low frequencies term, however, the above inequality fails. Now, remembering the 
definition of || • \b%^ ^ and the spectral properties of operator A_i , one has that 

||A_iVn||B^,, < c||A_iVn||Lco; 

at this point, Bernstein's inequality allows us to write that 

||A_iVn||sg^^^ < c ||vn||^2 . 

So putting together (fT7|) and the above inequality, we obtain 

(18) IIVHb^^^ < c (iivnii^. + ||An||^^-i 

First of all, let us see how to control ||An||os-i . Recall the third equation of (fTTI) : 
div (a Vn) = F with F := div {f - u ■ VVu). 
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Developing the left-hand side of this equation, we obtain 

F 

(19) An = -v(ioga) • vn + -• 

a 

Let us consider the first term of the right-hand side of the previous equation. 
If s > 1 then one may use that is an algebra and bound ||V(loga)||^s-i with ||Va||^s-i 

according to Proposition [3l we get 

||V(loga) • VnlLs-i < C||Va||RS-i||Vn||RS-i. 

_N_ _N_ 

Now, as Boo%o (see Corollary [2]) and B''^]. is an intermediate space between -Boo,oo and 

B^ j., standard interpolation inequalities (see e.g. [I], Chap. 2) ensure that 

(20) ||Vn||^s-i < C||Vn||^2||Vn||^^^^ for some ^ G]0,1[. 
Plugging this inequality in (jlSp and applying Young's inequality, we finally obtain 

|7 

I ds — 1 



(21) llVnllfi^^ < C (l + \\Va 





F 




) Iivn||i2 + 






a 


RS-l / 



where the exponent 7 depends only on the space dimension and on s. 

In the limit case s = r = 1, the space B^l is no more an algebra and we have to modify the 
above argument: we use the Bony decomposition ([6]) to write 

V(loga) • vn = rv(ioga)Vn + rvnV(loga) + i?(V(loga),Vn). 

To estimate first and second term, we can apply Propositions [2] and [3) we get 

(22) rv(ioga)Vn||so < c||v(ioga)||i«. iivniiso 

^ 00,1 00,1 

< C||Va||Loo||Vn||oo , 

00,1 

(23) ||rvnV(loga)||BO ^ < C ||Vn||L^ ||V(loga)||50 ^ 

< CllVnilioo llVallRO . 

CX3,1 

A similar inequality is no more true for the remainder term, though, nowever, one may use 

1 

that vn is in fact more regular: it belongs to B^ ^ for instance. Hence, using the embedding 

1. 

B^ ^ B'^ ^ and Proposition [2l we can write 

||i?(v(iogo),vn)||50^^ < C7||v(ioga)||L^ iivnii 1 , 

< c ||Va||L°° llvnii 1 . 

— II 11-^ II iIr^ 

Putting the above inequality together with ()22p and ()23p . and using that 1 we 
conclude that 

||v(ioga) • vniioo <c||Va||oo ||vn|| 1 . 

OC,l 00 .1 D ^ 

-°oo,l 

Now, using interpolation between Besov spaces, as done for proving (|2U|) . we get for some 
suitable 9 g]0, 1[, 

||v(ioga) • vnii^o < c||Va||so ||vnf ||vn||^2 . 

00,1 00,1 00,1 

Hence ||Vn||ni satisfies Inequality (I2ip for some convenient 7 > 0. 

00,1 

Next, let us bound the last term of (I19p . By virtue of Bony's decomposition ([6]), we have 

L = pF = TpF + Tfp + R{p,F) ■ 
a 

so from Proposition [2] we infer that 
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• \\TpF\\^s-i < Cp*\\F\\^s^,, 

• \\R{p,F)\\^s^^. < \\R{p,F)\\b^^^^ < CWpWb.^ JFW^.^-. < C\\p\\Bs^jF\\^s-y 

It is clear that ||div can be controlled by ^- For the second term of F we have to 

take advantage, once again, of Bony's decomposition ([6]) as follows: 

div (u-VVu) = Y,9tu'9,(Vuy = Y,(t9^u^9,Vu' + Tg^vu^d^u^ + diRiu^ , djVu')) , 



where in the last equality we have used also the fact that divT^u = 0. Now, for all i and j we 
have: 



s-1 

oo.r 



\diR{u^,djVu' 



Bt 



< C IIV-uIIlcx, ||V7^n| 

< C\\VVu\\lo^\\Vu\\jjs-i 



-1 < \\R{u^,djVu' 



IBS 



< C||u||rs IIVPuIIro 

II 11-^ run. II IIJ-'^-^ 



< CWuW 



\VVu\\ 



Because, by embedding, 
(24) 

we thus have 



iVPnllioo < C WWuW p,s-i , 



\\div{u-VVu)\\^s-i < C\\u\\bs^JVVu\\^s^^.. 

In order to bound Vu, let us decompose it into low and high frequencies as follows: 

Vu = A^iVu + (Id - A_i)Ptt. 

On the one hand, combining Bernstein's inequality and the fact that V is an orthogonal projector 
over yields 



\A^iWu\ 



< C \\u\ 



L2 



On the other hand, according to Remark [21 one may write that 

||(Id -A_i)Pn||B^„ <Ch||B^,, 

Therefore we get 



(25) 



from which it follows that 



(26) 



< C Hallos 



< C\\u\ 



-'oo.r 



Bi 



+ u 



B%o,r^L^ 



It remains us to control ||Vn||j;^2. Keeping in mind Lemma [21 from the third equation of 
System pT]) and Inequalities ([2l]) - ([25]) . we immediately get 

1,2 + \\u ■ Wu\\i2 
L2 + \\u\\i2 ||V7^m||loo 



||vn||i2 < 

< 



DENSITY-DEPENDENT INCOMPRESSIBLE EULER EQUATIONS 11 

Putting all these inequalities together, we finally obtain 

(27) ||Vn|L.(^.) < + jJ|Vll|g,,„nL^dr) 

(28) l|Vn||^j(^^^^) < C ((l + l|Va||;^(^^_.J ||Vn||^j(^.) + 

3.1.3. Final estimate. First of all, let us fix T > so small as to satisfy 

(29) """"^i,^ ii^ii^-,^^*) - 

a fact that is always possible because of the continuity of u with respect the time variable. 
Then, setting 

U{t) := ||n(t)||i2nBg^,, = h(t)||i2 + h(t)||Bs^„ 
Uoit) := ||no||L2nBg^^ + / ||/||L2nB^ . c^^- 







and combining estimates (dH), ([E]), ([I6D, ^ and ([28]), we get 

(30) U{t) < C {Uait) + ^ C/2(T)dr^ for aU t € [0,r], 

where the constant C depends only on s, A^, ||ao||_Bj^^, a* and a* . 

So, taking T small enough and changing once more the multiplying constant if needed, a 
standard bootstrap argument allows to show that 

u{t) < cuoit) yt e [o,r] . 

3.2. Existence of a solution to (jlip . We proceed in two steps: first we construct inductively 
a sequence of smooth global approximate solutions, defined as solutions of a linear system, and 
then we prove the convergence of this sequence to a solution of the nonlinear system (jllh with 
the required property. Recall that to simplify the presentation we have assumed that Tq = +oo 
and that we focus on the evolution for positive times. 

3.2.1. Construction of the sequence of approximate solutions. First, we smooth out the data 
(by convolution for instance) so as to get a sequence (oq , Uq , /")„gN such that Uq £ H°°, 
/" G C(IR+; -fr°°), Oq and its derivatives at any order are bounded and 

(31) a* < ao < a*, 
with in addition 

• < ao in B^,r: 

• — ^ in n B^ ,., 

. r^/in C(M+;L2)nLi(M+;i?^,,). 

In order to construct a sequence of smooth approximate solutions, we argue by induction. We 
first set a° = ag, it° = itg and Vn° = 0. 

Now, suppose we have already built a smooth approximate solution (a",?/", VII") over IR+ x 
with a" satisfying (I12p . In order to construct the (n + l)-th term of the sequence, we first 
define a"^^ to be the solution of the linear transport equation 

ata"+^ + n" • Va"+i = 

with initial datum a"+^|f=o = C'o'^^ ■ 



12 



R. DANCHIN AND F. FANELLI 



Given that is smooth, its flow is smooth too so that a""'"^(i, x) 



-.n+l 



((■0") ^{x)) , where 



■0" is the flow at time t. Note that ip^ is a smooth diff'eomorphism on the whole M^. From this 
fact, we gather that a"^^ is smooth and satisfies (jl2p . Furthermore, by virtue of Proposition U 



(32) 



-,n+l 



oo,r — II U ii-^oo,r 



exp C / ||n"||B^ dr 



Note that the reciprocal function p^'-^^ of a"+-^ satisfies x) = p^'^^ \{i>t) ^ 

together with (fT2]) and the equation 

Hence it also fulfills Inequality (p2|) up to a change of Cq"^^ in Pq^^. 

At this point, we define vT'^^ as the unique smooth solution of the transport equation 



n+l 



n+l 



-.n+l 



,n+l 



Since the right-hand side belongs to L\ 



div vl'dx^ I • 



n „,n+l 



dx. 



gf-x-^^+'i L"^) , from classical results for transport equation 
we get that n"^^ G C(M+;L^). Besides, as /)" = (a")^^ for all n, if we differentiate with respect 
to time the product \/ p^+^u""'"^ and take the scalar product with u^~^^ , we obtain 

1 

L2 ~ 

Observe that and u""*""*^ need not be divergence free; nevertheless one may control ||divtt"||Loo 
with ||ti"||_Bg^ ^ . So, from the previous equality, applying Gronwall's Lemma, it is easy to see 
that 

(33) ^ 



< 



L2 



+ c 

Finally, we have to define the approximate pressure n""*"^ . We have already proved that a"^^ 
satisfies the ellipticity hypothesis p^ : so we can consider the elliptic equation 

div (a"+i Vn"+i) = div (r+i - n"+^ • VPu"+^) . 

As f^~^^ and u""''^ are in C(M+; H°°), the classical theory for elliptic equations ensures that the 
above equation has a unique solution VII"'"'"^ in C(M+;-fr°°). In addition, going along the lines 
of the proof of ([27|) . we get 

Li{L2) + / II" llBi,.nL2 



(34) ||vn"+i||^,(^,) < C(||r+i|Li(^.) + ^ Il^x'^+MIL ...dr 



Of course, by embedding, we have Vn"+^ E C(M+; B^ ^.). Hence, arguing as for proving (f28]l . 
we get 

r-t 

(35) ||vn"+i| 



+ C 1 + llVa 



n+l ||7 



rn+li 



Note also that the norms of the approximate data that we use in ([32]) . (j33|) . (pH) and (j35]) 
may be bounded independently of n. Therefore, repeating the arguments leading to (pOj) and 
to Theorem 1 of [9], one may find some positive time T which may depend on UpollBgo^) 
IIuqIIb^ ^nL2 and ||/||li([o,t];B^ ^nL2) but is independent of n such that 

• (a")„eN is bounded in L°°([0, T]; S^^,), 

• (ti")„eN is bounded in L~([0, T]; S^'^^ n L^), 

• (Vn")„gN is bounded in L1([0, T]; J n L°°([0, T]; L^). 
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3.2.2. Convergence of the sequence. Let us observe that the function a"" := — satisfies 

J dta" = - u"-i • Va" 
\ a%=o = 0. 

Because m""^ G C{[0,T]; L"^) and Va"" G Cb{[0,T] x M^), it immediately follows that a" G 
C"^([0, r]; L^). Now we want to prove that the sequence (a", li", Vn"')„gN , built in this way, is a 
Cauchy sequence in C([0, r];L^). So let us define 



5" = , 





:= a"+i - 


a", 




:= - 


5" 




:= - 


P\ 




:= - 


n\ 




:= n"+i - 




<5/" 


:= - 





Let us emphasize that, by assumption and embedding, we have 

• Oq ao in C°'^, 

• Uq ^ uq in L^, 

. r ^ / in Ln[0,r];L2). 

This will be the key to our proof of convergence. 

Let us first focus on a". By construction, (5a" belongs to C^{[0,T]; L'^) and satisfies the 
equation 

atS" = -n" • V(E" - Su""^ ■ Va" - n" • V&[J 
from which, taking the scalar product in with S" , we obtain 

So, keeping in mind that 5z"(0) = and integrating with respect to the time variable one has 

(36) W\t)\\^2 < ||div^x"||^^ ||S"||^2 + ||Va"||i.o||&x"-i||i2 + ||^x"||i2 ||V&^||loo j dr. 

Equally easily, one can see that the following equality holds true: 

p"+i (at&i" + n" • V(5u") + Vffl"-! = (5/" - • Vn" - Vn"-^) ; 

taking the scalar product in with du"' , integrating by parts, remembering the first equation 
of (dl) at (n + l)-th step, we finally get 



L2 



< / lldivii'' 



L2 



dT 



+ 



|V^/"| 



+ 



L2 



L2 



+ 



lL2 



L2 



From (j36p . Gronwall's Lemma and (j3ip . we thus get for some constant C depending only on a* 
and a*, 



(5u 



n-l| 



IL2 



+ ||Vffl"-i||^2 + ||Vn"-i^2 II&oIIl- + \\u^L2\\VSa^\\Lo^ + ||5r||i2 ) dT, 



-rn—l 
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where we have set 

dr. 



:= ^ (||divn"||i^ + ||v^Vn"-i 



Of course, the uniform a priori estimates of the previous step ahow us to control the exponential 
term for all t E [0, T] by some constant Ct- 

Next, we have to deal with the term Viffl"""'^ . We notice that it satisfies the elliptic equation 
-div (a'^-^V-H"-^) = div (-(fa"-ivn" - u'^-^ • VVda''-^ - da"^'^ ■ VPu" + 5^-^) . 
Then applying the following algebraic identity 

(Mv{v-Vw) = div (tt; • Vw) + dvv {v divw) — div div 
to V = u"'^^ and w = T'du^~^ , and remembering that divVdu^~^ = 0, we get 

div (a"-^V(II"-i) = div (vSu''-^ divu^-^ - VSu''-^ ■ Vn"-i - (k""^ • VPn" 



Therefore, from Lemma[2]and the fact that ||'P||£(l2;L2) = Ij one immediately has the following 
inequality: 

(37) a,\\vm^-^\^, < ||fi"-i^2l|vn"||^^ + ||&^^||^^||vn"||^2 + 

+ \\6u^-'\\^, (||divn"-i||^^ + ||Vn"-i||^^ + ||VPu"||^^) . 
Due to a priori estimates, we finally obtain, for all t G [0, T] , 

||(<E",5a")(t)||^. < CT(||(5u^||L2+^*(||(&"-i,&i"-i)||^, 

+ ||Vffl"-i||^, + ll&^lbo,! + \\6n\L^) dr^ 

so, plugging the second inequality in the first one, we find out that for all t E [0, T], 

(38) \\{&^,Su^m\\L^<en + CT r||(S"-\&x"-l)||^, 

Jo 

with 

En := CT(^\\du^\\L2+ j\\\Sr-^L^ + WSPh^ + W&i'^-'Wl^ + ||V<5a^i||Loo) dty 
Now, we have 

y^En <00. 



From this and (I38p . it is easy to conclude that 

V sup (||5i"(t)||z.2 + ||&."(t)||z.2) <oo. 

n <e[0,T] 

In other words, {a'^)n<=N and {u'^)n<=N are Cauchy sequences in C([0,T];L^); therefore they 
converge to some functions a, u £ C([0,T];L^). In the same way, it is clear that (Vn")^gj^ 
converges to some VII E C([0,T];L^). 

Defining a := a + ao , it remains to show that a, u and VII are indeed solutions of the initial 
system. We already know that a, u and Vn E C([0,T];L2). In addition, 

• thanks to Fatou's property in Besov spaces, as (a")„gisj is bounded in L°°([0, T]; j.) , 
we obtain that a E L~([0, T]; ^) and satisfies i^; 
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• in the same way, u € L°°([0, T]; i?^ ^) because also {u'^)neN bounded in the same 
space; 

• finally, VII G L^([0, T]; i?^ ^) because the sequence (Vn"')^gpj is bounded in the same 
functional space. 

By interpolation we get that the sequences converge strongly to the solutions in every interme- 
diate space between C([0,T];L^) and C([0, T]; i?^ ^) , that is enough to pass to the limit in the 
equations satisfied by (a", it", VII"') . So, (a, u, VII) satisfies System (fTT]l . 

Finally, continuity properties of the solutions with respect to the time variable can be re- 
covered from the equations satisfied by them, using classical results for transport equations in 
Besov spaces (see Proposition H]). 

3.3. Uniqueness of the solution. Uniqueness of the solution to System ([1]) is a straightforward 
consequence of the following stability result, the proof of which can be found in [9j. 

Proposition 5. Let (/9i,ui,Vni) and (p2,U2,'VIl2) satisfy System ([T]) with external forces fi 
and f2, respectively. Suppose that pi and p2 both satisfy p^ . Assume also that: 

• 5p := p2 — pi and du := U2 — ui both belong to C^{[0,T]; L'^) , 

• 5/:=/2-/i GC([0,r];L2), 

• Vpi, Vui and VHi belong to L^([0, T]; L~) . 
Then for all t € [0, T] we have 



L2, 



< 



L2 



+ 



-""^'^ \\iVJ^2df)h. dr 



with 



A{t) :- 





Vpi 


+ 


VBi 


a 










Pl\/P2 



+ I|v^iIIl°° 



L°° 



Proof of uniqueness in TheoremUi Let us suppose that there exist two solutions (pi,ni,Vni) 
and (p2) '"2) Vn2) to System ([T]) corresponding to the same data and satisfying the hypotheses of 
Theorem [TJ Then, as one can easily verify, these solutions satisfy the assumptions of Proposition 
[5l For instance, that 6p G C^([0, T]; L"^) is an immediate consequence of the fact that, for i = 1,2, 
the velocity field m is in C{[0,T]; L'^) and Vpi is in C([0, T]; L°°) , so that dtpi £ C{[0,T]; L"^) . 
So, Proposition [5] implies that (pi, ui, VIIi) = (^2, ""2, VII2) . □ 

3.4. Proof of the continuation criterion. Now, we want to prove the continuation criterion 
for the solution to ([I]). We proceed in two steps. As usual, we will suppose Condition (C) to 
be satisfied with p = 00. The first step of the proof is given by the following lemma. 

Lemma 4. Let (p, n, VII) be a solution of System ^ on [0, r*[xM^ such tha^ 
. u€Ci[0,T*[;B'^.,)nC\[0,T*[;L^), 
• /9 G C([0,r*[;5^_J and satisfies (fT2]l . 

Suppose also that Condition (|4]) holds and that T* is finite. Then 

rT* 



sup 

tG[0,T* 



\U 



Bi. 



.nL2 + 



+ 



ivniifia dt < 00. 



Proof of Lemma [7} It is only a matter of repeating the a priori estimates of the previous section, 
but in a more accurate way. Note that a := 1/p satisfies the same hypothesis as p, so we will 
work without distinction with these two quantities, according to what is more convenient to us, 
and set g = p or a. Recall that 

dtq + u ■ Vq = 0. 

Hence, applying operator Aj yields 

dtAjq + u ■ VAjq = [n • V, Aj]q 



with the usual convention that continuity in time is weak if r — 00. 
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whence, for all t £ [0,T*[, 



(39) 



2^^||A,-(7(t)||ioo < 2^^||A,go||L- + / 2^1| [n • V, A,]g||Loo dr. 

^0 



Now, Lemma 2.100 in [1] ensures that 



(2^'^|| [u ■ V, Aj]q\\L^)^ < C{\\Vu\\l^ h\\Bi^^^ + IIV^zIIl- I|V^x||^^-i). 



dr 



Hence, performing an £^ summation in (j39p . we get 

(40) \\q{t)\\Bs^, < Ikobj^, + C [ {\\Vu\\L^\\q\\Bs^^ + \\Vq\\L^\\u\\Bs^^)dT. 

Jo 

As regards the velocity field, we have according to (jl4|) . 

Mt)\\L^<c(^\\uo\\L2 + j\\f 

and the last part of Proposition H] guarantees that 
\\u{t)\\Bi^^^ < exp j ||Vu||loo dr^ 

Bounding the last term according to Corollary [3l we thus get 

\\u{t)\\Bi^^^ < exp j \\Vu\\l^ dr^ 

X (ll^olb^,. + l\-Clo\\^uU^dr' ^ii^ii^^^^ ^^.ji^nbg^,, + ||Va||^^-.||Vn||i^) dr 
As regards the pressure term, we have 

||vn|U2 < c{\\f\\L2 + \\u\\l2 \\Vu\\l^) 



Rs + llaVnll Ra 



dr 



ivni 



oo,r 



< c ||vn||r2 + iiva- vniLs^i + 

III iJr^nr 



1 



div {f — u - Vu] 



Note that Bony's decomposition combined with the fact that div u = ensures that 

||div {u ■ Vu)\\^s-i < C||Vti||Loo||n||Bg^_^. 
In addition, under the assumption that s > 1 , Corollary [3] implies that 
(41) ||Va • Vn||^.-i < C(||Va||Loo||Vn||gs-i + ||Va||^.-i||Vn||ioo). 



So finally 



Ivni 



<c{ ||vn||i2 + ||Va||Loo||vn||^s_-i + ||Va||Rs-i||vn 



IBS 



\L°° 



Rs + II Vnll roo ||n|| Rs 



Putting together all these estimates and applying Gronwall's Lemma, we obtain if s > 1, 
||Va||^j-i + ||n(t)bg^^^nL2 < «^p(^_{ ||(Va, V^x, Vn)||i^(ir^ (||Vao||B^-i + 



l-wollij^^nL^ + ||/||ij^^nL2 + / l|Va| 

10 



ivni 



dr 



where the constant C depends only on s, a^,, a* and N. 
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Now, the equation for Va and Gronwall inequahty immediately ensure that 
(42) ||Va(t)||Loo < ||Vao||L°° exp \\Vu\\l^ dr^ , 

which, thanks to Hypothesis (jl|) implies that Va is bounded in time with values in L°° . 

Moreover, by hypothesis VU G L^{[0,T*[; B^^l) and Vu £ L^{[0,T*[; L°°); at this point, 
keeping in mind the embedding B^}. ^ L°° , the previous inequality gives us the thesis of the 
lemma in the case s > 1. 

In the endpoint case s = r = 1, Inequality (I4ip fails. In order to complete the proof of 
the lemma, we will have to take advantage of the vorticity equation associated to ([I]). This is 
postponed to the next section. □ 

The second lemma, which will enable us to complete the proof of Theorem [2] reads: 

Lemma 5. Let (p, u, VII) be the solution of System ([1]) such tha^ 
. pGC([0,T*[:B^^,,) and (HI]); 
. nGC([0,r*[:i?^,,)nCi([0,r*[;L2); 
. vn e C([0, T* [; L2) n lH[o, T*[; i?^,,). 
Moreover, suppose that 

Ihlli- {iJ|„,,nL2) + l|Va||^co (5^-1) < oo. 

Then (p, u, VII) can be continued beyond time T* into a solution of ([1]) with the same regu- 
larity. 

Proof of Lemma\^ From the proof of Theorem[T]we know that there exists a time e, depending 
only on p*, s, \\u\\ [b^ ^nL'^) ^ II^o|Il°o (b""^) norm of the data such that, for 

all T < T , Euler system with data {p{T),u{T), f(T + •)) has a unique solution until time e. 

Now, taking for example T = T — e/2, we thus obtain a solution, which is the continuation 
of the initial one, (p, n, VII) , until time T + e/2. □ 

Let us complete the proof of Theorem [2j The first part is a straightforward consequence 
of these two lemmas. Indeed: Lemma |4] ensures that (B^ ^nL^) and ||Va|| are 

finite. As for the last claim (the Beale-Kato-Majda type continuation criterion), it is a classical 
consequence of the well-known logarithmic interpolation inequality (see e.g. |1]) 



|Vii||L=o < C \\u\\i2 + ||ri||Loo log e + 



ns-l 



|r2||/^oo / / 

So Theorem [2] is now completely proved, up to the endpoint case s = r = 1. □ 

4. The vorticity equation and applications 

This section is devoted to the proof of the blow-up criterion in the endpoint case s = r = 1, 
and of Theorem [3l Both results rely on the vorticity equation associated to System ([T|) . As done 
in Section [3l we shall restrict ourselves to the evolution for positive times and make the usual 
convention as regards time continuity, if r < cxd. 

4.1. On the vorticity. As in all this section the vorticity will play a fundamental role, let us 
spend some words about it. Given a vector-field u, we set Vu its Jacobian matrix and *Vti the 
transposed matrix of Vu. We define the vorticity associated to u by 

n := Vu - *Vu. 

Recall that, in dimension N = 2, Q, can be identified with the scalar function oj = diu^ — 
d2U^ , while for = 3 with the vector-field u = V x u. 



^with the usual convention that continuity in time is weak if r — oo. 
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It is obvious that, for all q G [l,oo], if Vu £ , then also Q £ L'^ . Conversely, if u is 



divergence- free then for all 1 < i < we have An' = X^jLi dj^ij-, and so, formally, 

N 

Vn' = -V (-A)-^ ^ dj Clij . 

i=i 

As the symbol of the operator —di{—A)^^dj is a{^) = CiCj/lCP; the classical Calderon- 
Zygmund Theorem ensures tha10 for all q G ]1, oo[ if G L'^ , then S/u G L'^ and 

(43) ||Vu||l, < C||J1||l,. 

The above relation also implies that 

u = A_iu - (Id - A_i)(-A)-^ ^i^v- 



Hence combining Bernstein's inequality and Proposition [H we gather that 
(44) Iklloi < C(||n||L,j + ||0|| oo ) for ah gG[l,ool. 

oo,l ^ 00,1'^ 



From now on, let us assume that VL is the vorticity associated to some solution (p, u, VII) of 
([1]), defined on [0,r] x M^. From the velocity equation, we gather that Q. satisfies the following 
transport-like equation: 



(45) 



+ u • vj^ + • Vn + *vu • + V ( i ) A vn = F 



where Fij := djP — dif^ and, for two vector fields v and w, we have set v A w to be the 
skew-symmetric matrix with components 

{v A w)^j = v^w^ — v'^w^ . 

Using classical L'^ estimates for transport equations and taking advantage of Gronwall's 
Lemma, from (145p we immediately get 



(46) \\^{t)\\L<i < '^^pi^J l|Vii||L°°dr 











^Vp A vn 











Let us notice that, in the case of space dimension N = 2, equation (I45p becomes 

1^ 



dtuj + u-Vu + V i^-j A vn = F , 
so that one obtains the same estimate as before, but without the exponential growth: 



|^^(t)||L'^ < II^^(o)||l'^ + 











4tV/9A vn 


J 




p 





dr. 



Therefore, the two-dimensional case is in a certain sense better. We shall take advantage of that 
in Section [51 As concerns the results of this section, the proof will not depend on the dimension. 
So for the time being we assume that the dimension is any integer greater than or equal to 2. 



'This time the extreme values are not included. 
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4.2. Proof of Theorem [2] in the limit case s = r = 1. We just have to modify the proof 
of Lemma m From the vorticity equation (j45p and Proposition U] (recall that divu = 0), we 
readily get 

ft 



(47) ||O(t)||50^^ < exp^^C^ llVullLoodr^ 



X ||17o||bo + / llFll^o dr + / (||VaAVn||^o + • Vn + 'Vu • f^ll^o jdr). 
We claim that 

(48) ||Va A Vnll^o^^ < C(||Va||Loo||Vn||50^ ^ + || VHUloo ||a||Bi^ J , 

(49) Vu + *Vu-il||oo < C||Vn||ioo||tt||gi . 

oo , 1 oo ,1 

Both inequalities rely on Sony's decomposition ^ and algebraic cancellations. Indeed, we 
observe that 

diadjR - djadiU = Tg^adjH - To^adiU + TQ^ndi.a - Tg^ndja + d,R{a, djii) - djR{a, diU). 

Applying Proposition [2] thus yields ()48p . 
Next, we notice that, as divn = 0, 

{n-Vu + ^Vu-n)^. = ^{diu'^dku^ -dju'^dku'), 

k 



Therefore, 

{n-vu + *vti • n) 



k ^ 

Hence Proposition [2] implies (|49p . 

It is now easy to complete the proof of Lemma S] in the limit case. Indeed, plugging ()48p and 
9|) in (jI7|) . using the energy inequality (fH|) and Inequality (jH]) with q = 2, we easily get 



||n(t)||gi ^ p|£2 < C exp ^(7 J \\\7u\\loo dr 

x(\\uo\\b^ ,nL2+ / ll/llfii ,nL^dT+ [ (||Va||L.c ||Vn||^o + || VHUloo ||a||Bi ) dr] . 
Hence, adding up Inequality (j^Oj) and using Gronwall's inequality, we end up with 
Xit) < Cexp(^C l\\iVu,Va,VU)\\Lo^ dr^ (^X(O) + ^'(||/||bi^ ^nL^ + II VaU^oo || VnH^o^ J dr^ 
with X{t) := \\a{t)\\j,i +\\u{t)\\Bi nL2. 

oo , 1 oo,l 

Taking advantage of (|42p completes the proof of Lemma H] in the limit case. □ 

4.3. Proof of Theorem [3l We first prove a priori estimates, and then we will get from them 
existence and uniqueness of the solution. In fact, it will turn out to be possible to apply Theorem 
[1] after performing a suitable cut-off on the initial velocity field and thus to work directly on 
System ([1]) , without passing through the equivalence with (jlip or with a sequence of approximate 
linear systems. 
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4.3.1. A priori estimates. As in the previous section, remembering also Remark [3l the fohowing 
estimates hold true: 

(50) l|Vp(t)IU^-i < ||VpolU-i exp(c^Vb^,.rf^) 

(51) hWllsg^,,, < exp^C^ \\u\\Bs^,.dT^ ■ (^\\uo\\b^^^, + 

^0 / 

Moreover, from the transport equation satisfied by the velocity field, we easily gather that 

dr. 



f 


vn 




Jo 


p 


L4 



U[t)\\Li < \\uo\\l^ + 



Therefore, using interpolation in Lebesgue spaces and embedding (see Corollary [2]) , 

||n(t)||i4 < ||uo||l4 + — / ||Vn|||o. llVnlll^dr 

P* Jo 

C /"* 

(52) < ||wo||l4 + — / \\VU\\Bs^ nL^dT. 

P* Jo 

In order to bound the vorticity in , one may use that 



1 



Vp A vn 



L4 



< 



< 



1 



{p*r 
1 

c 



||Vp A vn||^4 
IIVpIIl- ||vn||^4 



From this and ()46p . we thus get 

(53) \\n{t)\\L4 < exp^^ \\Vu\\^s-idT^ 



X l|f^o||L4 + 



C 



{P* 



Now, in order to close the estimates, we need to control the pressure term. Its Besov norm 
can be bounded as in Section [3l up to a change of \\u\\i2 into ||tf||L4; indeed it is clear that 
in Inequality (j25p the norm of u may be replaced by any norm with q < oo. As a 
consequence, combining the (modified) inequality (j26p and ()2ip yields 
(54) 

In order to bound the norm of VII, we take the divergence of the second equation of 
System ([T]). We obtain 



-div 



V P 



div (n • Vu) , 



from which, applying elliptic estimates of Lemma [2] and 
(55) ||Vu||i4 < C||f]||^4, 

we get 



(56) 



1 

p* 



VIl\\i2 < \\u-'Vu\\r,2 < l|u|| r4 II Vn|| r4 < C||u 



|L2 



L4 



|L4 



\l* II^IIl4 
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We claim that putting together inequalities §0^, 1^, 1^, §3^ and §6^) enables 

us to close the estimates on some nontrivial time interval [0, T] depending only on the norm of 
the data. 

In effect, assuming that T has been chosen so that Inequality ()29p is satisfied, we get from 
the above inequalities 

\\u 



l,^, < 2||'Uolbg,,,, +C'o||Vn||ii(B^ J, 



\\u{t)\\L4 < ||no||L4 + Co||Vn||ri/Bs ) + C'o/ II^^IIl^II^^IIlic?-^, 

Wmh^ < 2||no||i4 + Co||vn||ii(5^_^), 

where the constant Cq depends on s, p*, /?*, N and ||/Oo||ri • 

oo,l 

Therefore, applying Gronwall lemma and assuming that T has been chosen so that (in addition 
to ([29])) we have 

fT 

||u||;4/i,4 dr < c 





where c is a small enough constant depending only on Cq, it is easy to close the estimates. 

Remark 4. Exhibiting an L? estimate for VII even though u is not in is the key to the 
proof. This has been obtained in (|56p . Note however that we have some freedom there. In 
fact, one may rather assume that uq G and Vuq € L'^ , with p and q in ]2,cxd[ such that 
1/p + 1/q > 1/2 and get a statement similar to that of Theoreml^ under these two assumptions. 
The details are left to the reader. 

4.3.2. Existence of a solution. We want to take advantage of the existence theory provided by 
Theorem [TJ However, under the assumptions of Theorem [3l the initial velocity does not belong 
to L2. To overcome this, we shall introduce a sequence of truncated initial velocities. Then 
Theorem [1] will enable us to solve System ([1]) with these modified data and the previous part 
will provide uniform estimates in the right functional spaces on a small enough (fixed) time 
interval. Finally, convergence will be proved by an energy method similar to that we used for 
Theorem [TJ 



First step: construction of the sequence of approximate solutions. Take any $ S Cq°(R^ ) with 
$ = 1 on a neighborhood of the origin, and set ^nix) = ^(x/n). Then let us define Uq := uq 
for all n E N. 

Given that Uq is continuous and compactly supported, it obviously belongs to L^. Of course, 
we still have Uq £ ^ n VT^'^ H , so we fall back into hypothesis of Theorem [TJ From 
it, we get the existence of some time Tn and of a solution (p",n",Vn") to ([T|) with data 
{po,u^,0) such that G C([0,r„];B^^,), n" G Ci([0, r„]; L^) n C([0, r„]; S^^,) and VH" G 
C([0, Tn]; L"^) n -^^"'^([0, Tn]; B^ j.) . From (f55]l . the vorticity equation and the velocity equation, it 
is easy to see that, in addition, it" G C([0,T„]; W^'^). 

Finally, as the norm of Uq in W^''^ H B^ ^ may be bounded independently of n, the a priori 
estimates that have been performed in the previous paragraph ensure that one may find some 
positive lower bound T for T„ such that (/o", u", VII") satisfies bounds independent of n on 
[0, T] in the desired functional spaces. 

Second step: convergence of the sequence. As done in the previous section, we define = 
p"' — po, and then 



,5u" 
dl" 
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Resorting to the same type of computations as in the previous section (it is actually easier 
as, now, divti" = for all n), we can prove that {p^,u^,V'n."')^^^ is a Cauchy sequence in 
C([0,r];L2). Hence it converges to some {p, u, VII) which belongs to the same space. 

Now, defining p := pQ + p, bearing in mind the uniform estimates of the previous step, and 
using the Fatou property, we easily conclude that 

. p G L~([0,r];S^ ,,) and p. < p < p* ; 
. n G ( [0, T] ; i?^,,) n ( [0, T] ; . 

. Vn e Li([0,T];S^,,) nL°°([0,r];L2). 

Finally, by interpolation we can pass to the limit in the equations at step n, so we get that 
u, vn) satisfies ([T]), while continuity in time follows from Proposition HI □ 

5. Remarks on the lifespan of the solution 

In this section, we exhibit lower bounds for the lifespan of the solution to the density dependent 
incompressible Euler equations. We first establish that, like in the homogeneous case, in any 
dimension, if the initial velocity is of order e then the lifespan is at least of order even if 
the density is not a small perturbation of a positive real number. Next we focus on the two- 
dimensional case: we show in the second part of this section, that for small perturbations of 
a constant density state, the lifespan tends to be very large. Therefore, for nonhomogeneous 
incompressible fluids too, the two-dimensional case is somewhat nicer than the general situation. 

5.1. The general case. Let pQ, uq and / satisfy the assumptions of Theorem [1] or [3j Denote 

uq{x) := e~^uo{x) and f{t,x):=e^'^f{e^^t,x). 
It is clear that if we set 

(/), u, vn) (t , x) = {p, eu, e^Vn) {et, x) 

then (p, 5, vn) is a solution to ([1]) on [T^,T*] with data {pq,uq, f) if and only if (/>, n, Vn) is 
a solution to ([T|) on [e^-'^T^,, e^^T*] with data {p(j,U(j,f). 

Hence, putting together the results of the previous section, we can conclude to the following 
statement. 

Theorem 4. Let {pq,uq) satisfy the assumptions of Theorem{l\ orl^ and / = 0. There exists a 
positive time T* depending only on s, N, p^,, ||/Oo||ro and ||5o||ro such that for any e > 

oo,l oo,l 

the upper bound T* of the maximal interval of existence for the solution to ^ with initial data 
{Po,£Uq) satisfies 

T* > £-^T*. 

A similar result holds for the lower bound of the maximal interval of existence. 

5.2. The two-dimensional case. Recall that for the homogeneous equations, any solution 
corresponding to suitably smooth data is global, a fact which relies on the conservation of the 
vorticity by the flow. Now, in our case, the vorticity equation reads (if / = 0) 

(57) dtUJ + u-Vuj + VbAVU = 

with 6 := 1/p - 1 and V6 A VH := dibd2U - d2bdiU. 

Owing to the new term involving the pressure and the nonhomogeneity, it is not clear at all 
that global existence still holds. Nevertheless, we are going to prove that the lifespan may be 
very large if the nonhomogeneity is small. 

To simplify the presentation, we focus on the case where po G B^^ li^"^) and uq G 5^ i(I^^) 
(note that Corollary [1] ensures that this is not restrictive) and assume, in addition, that uq G 
H^{M.'^) (this lower order assumption may be somewhat relaxed too). 

We aim at proving the following result. 
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Theorem 5. Under the above assumptions there exists a constant c such that if bo := — 1 
satisfies 

(58) ll^ollsi <c 



oc ,1 



then the lifespan of the solution to the two-dimensional density dependent incompressible Euler 
equations with initial data {pq,uq) and no source term is bounded from below by 

71 — n— log ( 1 + loe 



FollHinBi , V \\"o\\b 



Proof. Let ]T^,T*[ denote the maximal interval of existence of the solution (p, u, VIl) corre- 
sponding to {po,uo)- To simpUfy the presentation, we focus on the evolution for positive times. 

The key to the proof relies on the fact that in the two-dimensional case, the vorticity equation 
satisfies ([57|) . Now, it turns out that, as discovered by M. Vishik in [15] and by T. Hmidi and S. 
Keraani in [12j, the norms in Besov spaces with null regularity index of solutions to transport 
equations satisfy better estimates, namely in our case 

\\oj{t)\\BO^_^ <c(^\ujo\\bo^^+ J A Vnll^o^ ^ dr^ ^1 + ^ \\Vu\\LoodT^ 

whereas, according to Proposition HI the last term has to be replaced with exp^ J \\Vu\\l°° dr^ 

for nonzero regularity exponents. 

Therefore, using Inequality (jl8|) . we get 

(59) \\uj{t)\\BO^^^<c(^\\uJo\\BO^^^+ \\b\\BlJVU\\BO^^^dT^(^l + WVuU^dT^ 

Of course, a basic energy argument leads to 

(60) \\uj{t)\\L2 <\\ujo\\l^ + [ ||V6||Loo||Vn||i2dr 

Jo 

and it is well-known that for two-dimensional divergence-free vector fields, we have 

Therefore putting together Inequalities ([59]) and ([60|) and bearing in mind Inequality (ji^ and 
the energy inequality for u, we get 

(61) X{t) <c(^Xo + B\\VU\\so^^^L^dT^(^l + XdT^ 
with 

Xit):=Mt)\\H^nBl^, and B{t) := \\b{t)\\s.^^^. 
Bounding B is easy given that 

dtb + u-Vb = 0. 

Indeed, Inequality d?]) ensures that 

\\m\\B^ , < II&oIIbi ^expfc /V^IIbo .dr). 

00,1 00,1 y 00,1 J 

Therefore, 

(62) B{t) < Bo exp ^ ■ 

Bounding the pressure term in B'^ L'^ is our next task. For that, recall that, as 

div (aVn) = -div {u ■ Vu), 
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Lemma [2] guarantees that 

(63) a*||Vn||2,2 < ||u||j^2 II V-uIIloo . 

Next, differentiating once the pressure equation and applying again an energy method yields 

(64) a*||V^n||i2 < ||Vu||i2||Vu||Loo + ||Va||Loo || Vn||i2 . 
Therefore, combining (j63p and ([64]) and using elementary embedding, we get 

(65) l|Vn||^i < C\\a\\Bi^^X'^. 

Note that ||a|| ri and A:= \ + B are of the same order. This will be important in the sequel. 

oo,l 

In order to bound the pressure term in ^ we shall use the following classical logarithmic 
interpolation inequality (see e.g. [1], Chap. 2): 



ivni 



(66) llVnllso <C||Vn||^ilog e + 

°°' V llv^J-ll/fi 

In order to estimate ||Vn|| ni , we use the identity 

CXD,1 

Vn = A_iVn + A{D)diy {u ■ Vn) + A{D)diy (6Vn). 
with A{D) := (-A)-iV(Id - A_i). 

On the one hand, combining Bony's decomposition with the fact that div {u ■ Vn) = Vn : Vn, 
it is easy to show that 

||div (n • Vn)|| no < ||n||Ri . 
On the other hand. Proposition [2] guarantees that 

llfevniisi ^ <c(||6||Loo||vn||Bi ^ + II vniiioo 11611^1 

CXD , 1 \ CXD , 1 OO , 1 



So using the fact that A{D) (resp. A{D)div ) is a multiplier of degree —1 (resp. 0) away 
from the origin, we get from Proposition [H 



livnii^i < c ||vn||i2 + iinirRi + iivniiLco 11611^1 + ||6||Lco||vn||5i 

oo,l \ OO 1 00,1 OO 

Note that ||6(t)||ioo is time independent and that B^^ ^ L°° . Hence, under assumption (f58l) 
with c small enough, the last term may be absorbed by the left-hand side. As regards the last 
but one term, we use the following interpolation inequality: 

||vn||Lco < C7||vn|||2||vn|||i . 

oo,l 

Combining with Young's inequality, we thus conclude that, under assumption ([58]) . we have 

livniigi ^<c(\\u\\l, +{i + \\h\\l, )||vn||^2 
00,1 \ -"00,1 -"00,1 



Bounding the last term according to ([63|) . we thus end up with 

livniioi <CA^x'^. 

Inserting this inequality in (|66p and using also (j65p . one may now conclude that 

(67) llVnbo^^^nL^ < CAX^ log(e + B). 

It is now time to insert Inequalities (I62p and (I67p in (I6ip : we get 

(68) X{t)<c(^Xo + BoAolog{e + Bo) ^ X^ dr^ (^1 + X dr^ . 

Let To denote the supremum of times t £ [0, T* [ so that 

i-t 

(69) BoAolog{e + Bo) / X^ dr < Xq. 

Jo 
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From (j68p and Gronwall's Lemma, we gather that 

X{t) < 2C7Xoe2^*^« for all t G [0,ro[. 
Note that this inequality implies that for all t G [0, Tq [, we have 

Therefore, using ()69p and a bootstrap argument (based on the continuation theorems that we 
proved in the previous sections), it is easy to show that Tq is greater than any time t such that 



^0^0 log(e + So) (^e4^*^« - 1^ e^v{c{e 



2CtXo 



1 < 1. 



Taking the logarithm and using that logy < y — 1 for y > 0, we see that if Bq is small enough 
(an assumption which implies in particular that j4o log(e + Bq) is of order 1) then the above 
inequality is satisfied whenever 



C + 2 \2CBo 

This completes the proof of the lower bound for T* . □ 

Remark 5. If ojq has more regularity (say ujq G for some r G (0, 1) ) then one may first write 
an estimate for ||a;||j;^oo and next use the classical logarithmic inequality for hounding ||Vn||/,oo 
in terms of \\(jj\\loo and ||cL!||cr. The proof is longer, requires more regularity and, at the same 
time, the lower hound for the lifespan does not improve. 
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